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The interpretation of the k dependent spectral functions of the one-dimensional, infinite U Hubbard 
model obtained by using the factorized wave-function of Ogata and Shiba is revisited. The well 
defined feature which appears in addition to low energy features typical of Luttinger liquids, and 
which, close to the Fermi energy, can be interpreted as the shadow band resulting from 2&f spin 
fluctuations, is further investigated. A calculation of the self-energy shows that, not too close to 
the Fermi energy, this feature corresponds to a band, i.e. to a solution of the Dyson equation 
w-e(fe) -ReE(k,w) = 0. 
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In the study of systems of interacting electrons, the 
dynamical correlation functions play a central role be- 
cause they correspond to the response functions that can 
be measured in several experiments (infrared reflectivity, 
inelastic neutron scattering, nuclear magnetic resonance, 
etc..) The most basic of these correlation functions is 
the time-ordered Green's function, which can be defined 
by its spectral representation according to 
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The spectral functions A{k,u) and B{k,u) are defined 
by 

A(k,u J )=J2\(f,N+l\c{jO,N)\ 2 5(u;-E? +1 + E») 

f,a 

B(k,u;) =J2\(f,N-l\c k j0 7 N)\ 2 Siui-Elf + E?- 1 ) 

f,<r 

in the standard notation, and in principle they can be 
measured in angular-resolved inverse photoemission and 
photoemission experiments, respectively. The effect of 
electron-electron interactions on the spectral functions 
of the three-dimensional Coulomb gas has been investi- 
gated in much detail several years ago H . With respect 
to the simple 5(u> — £&) structure of the non- interacting 
case, the main differences are: i) A shift of the energy 
of the quasi-particle band; ii) A broadening of the quasi- 
particle peak; iii) A renormalization of the weight of the 
quasi-particle peak, compensated by the appearance of 
an incoherent background; iv) The presence of another 
band - a plasmon band - due to the long-range nature of 
the Coulomb potential. For a short-range repulsion, such 
as the on-site interaction of the Hubbard model, there is 
no plasmon band, and the spectral function is expected 
to have only one well-defined feature, the quasi-particle 
band, on top of an incoherent background. The validity 
of this simple picture in lower dimensional systems is cur- 
rently under intense discussion. In two dimensions (2D), 



there are at least two important issues: 1. the actual ex- 
istence of a quasiparticle peak, and 2. the generation by 
antifcrromagnetic short-range fluctuations close to half- 
filling of additional well-defined features in the spectral 
functions, called generically "shadow bands" 0. The ab- 
sence of exact results in 2D for correlated models makes 
the interpretation of the experimental results quite diffi- 
cult. For high-T c cuprates, these issues are still contro- 
versial ||. 

In ID, a number of analytical approaches are available. 
For instance, models having both charge and spin low- 
lying excitations belong to the universality class of the 
Luttinger liquid, and bosonization gives an accurate de- 
scription of their low-energy properties. The consequence 
for the spectral functions is that there is no quasi-particle 
eak but two divergences due to spin-charge separation 
. A lot of insight has also been obtained by the Bcthc 
ansatz solution of integrable models, the prototype being 
the Hubbard model defined by: 
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where c] „ and c„ „ arc electron creation and annihilation 



operators, and rii t(7 — c[ a c l a 



is the density operator. The 
Bethe ansatz wave-functions are so complicated that it 
is not possible to use them to calculate correlation func- 
tions. However, Ogata and Shiba showed a few years 
ago that, in the limit U — » +oo, the ground-state wave 
function can be written as a product of a spinless fermion 
wave-function \ip) and a squeezed spin wave-function \x), 
and that this wave function allows a very precise evalu- 
ation of the momentum distribution function H. More 
recently, Penc et al [DQ] used the same representation 
for the excited states to calculate the spectral functions 
A(k,u) and B(k,ui). The results for n = 1/4 are shown 
in Fig. 1. 

There are several interesting features to notice. In the 
low energy region near kp we can identify three struc- 
tures. For k < kp there are maxima at u = u c (k — kp) 



and uj ~ and a lot of spectral weight between them. 
There is also a small weight appearing on the other side 
of the Fermi energy for lu > —u c (k — kp). If we re- 
member that the spin velocity u s vanishes for the in- 
finite U Hubbard model, all these features are qualita- 
tively consistent with the Luttinger liquid calculations 
H . The dispersion of the charge part is exactly given by 
E(k) = -2tcos(\k\ + k F )- 




FIG. 1. Spectral functions of the U —* +00 Hubbard model 
for L = 228 sites and N — 114 electrons (n = 1/4) with Fermi 
momentum Hf = 7r/4. 

However, the Luttinger liquid picture does not exhaust 
the features of the spectral function of Fig. 1. For larger 
energies, or away from k F , there is a well defined band- 
like structure with considerable spectral weight and a 
dispersion given by E(k) = — 2icos(— \k\ + kp). This 
band crosses the Fermi energy at 3kp and has been inter- 
preted as a shadow-band originating from the diverging 
2k F spin fluctuations. The presence of a shadow-band in 
a ID lattice model has been first reported on the basis 
of exact diagonalizations of finite clusters by Haas and 
Dagotto g . More recently, a systematic study of various 
extensions of the infinite U Hubbard model to finite re- 
pulsions has shown that these features are robust down 
to reasonable values of the parameters [pi . 



While there is little doubt that this physical picture 
is essentially correct, the interpretation should be made 
more precise in two respects. First of all, the concept of 
a shadow band is no longer well defined below the energy 
where the charge part of the main band crosses k = 0. 
In fact, if we think in terms of the whole Brillouin zone 
(— 7r < k < 7r), what we call a shadow-band for k > is 
continuously connected to the main band for k < 0. So 
what we have called a shadow band is probably better 
thought of as a band per se which, close to the Fermi 
level, can be interpreted as a shadow band. 

Second, the term band has a precise meaning in the 
context of perturbation theory. A feature of the spec- 
tral function is called a band if it corresponds to a solu- 
tion of the Dyson equation u> — e(k) — Re£(k, lu) = 0, 
where e(k) = —2tcosk is the non-interacting disper- 
sion while the self-energy E(fc,w) is implicitly defined by 
G(k,co) — l/(ui — e(k) — T,(k,uj)) and can be obtained 
from G(k,uu) as T,(k,uj) = u — e(k) — G~ l (k,oj). For 
ID systems, the perturbation expansion diverges notori- 
ously, and, even for small U/t, T,(k,uj) cannot be calcu- 
lated perturbatively. However, having obtained A(k,ui) 
and B(k,tu) in an essentially exact way, we can use Eq. 
(1) to calculate the time ordered Green's function, and 
from it the self-energy The results for the real part and 
the imaginary part of the self-energy for k = 7r/3 and 
k = 2ir/3 are given in Fig. 2. 

To discuss the relationship between the spectral func- 
tions and the self-energy, it is useful to write A(k, uj) and 
B(k, uj) in terms of £(fc, lu) as 
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Let us start with k = 2ir /3. As for a non-interacting 
system, the Dyson equation has only one solution cor- 
responding to the intersection of the straight line with 
the real part of E(fc,w). This corresponds to the charge 
part of the main excitation band. The spin part, which 
is pinned at the Fermi energy in the present case, cor- 
responds to a vertical slope at lu = + in the imagi- 
nary part of the self-energy, as well as in the spectral 
function A(k, to). The "shadow" band, which appears 
at negative energies for that wave-vector, does not cor- 
respond to a solution of the Dyson equation, but to a 
maximum of the imaginary part that occurs in a region 
where uj — e(k) — Re£(k, lu) is large, so that the denomi- 
nator of A(k,uu) is dominated by (uu — e(k) — ReS(k, uj)) 2 . 

For k — 2tt/3, there are three solutions to the Dyson 
equation. As in the previous case, the solution for posi- 
tive uj corresponds to the charge part of the main band, 
and the spin part also gives rise to a singularity at 
lu = + . For negative frequencies, the situation is very 
similar to that of the plasmon band of the 3D Coulomb 
gas pF. Decreasing uj starting from the Fermi energy, 
the first solution of the Dyson equation coincides with 



a maximum of ImS(k, u>) and does not give rise to any 
feature in the spectral function. On the contrary, the 
other solution occurs in a region where Im£(k, u>) is not 
particularly big. This solution gives rise to a feature in 
the spectral function that corresponds to what we have 
called the "shadow" band. 
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case however it can't really be interpreted as a "shadow" 
band because it cannot be connected to the main band 
by a vector 2fci?. Close enough to the Fermi energy, this 
feature can really be thought of as a shadow of the main 
band because it really follows the charge part of the main 
band at a distance 2fcp apart. However, this feature does 
not correspond any more to a solution of the Dyson equa- 
tion, and strictly speaking, it should not be called a band. 

To summarize, this work is the first attempt at us- 
ing the self-energy to interpret the spectral functions of 
the Hubbard model that have been obtained in an es- 
sentially exact way in the U — ► +oo limit. Extracting 
the self-energy from the spectral functions numerically 
has already allowed us to draw an interesting conclusion 
concerning the origin of the feature that was previously 
interpreted as a shadow band, namely that this feature 
only corresponds to a solution of the Dyson equation 
far enough from the Fermi level. The next step, a de- 
tailed analysis of the behaviour of the self-energy close to 
(lu = 0,k = kp, 3kp), is in progress. It might help clarify 
the still controversial problem of the spectral function in 
2D in the presence of almost perfect nesting. 

We thank H. Shiba and A.-M. Tremblay for very useful 
discussions. 
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FIG. 2. Frequency dependence of the spectral functions 
and of the real and imaginary parts of the self-energy of the 
quarter-filled, U — * +00 Hubbard model for selected values of 
the momentum. The straight line corresponds to uo — e(k). a) 
k = tt/3; b) k = 2tt/3. 

Coming back to the second point raised above, we see 
that the feature that crosses the Fermi energy at 3fcp 
corresponds to a band, i.e. to a solution of the Dyson 
equation, only far enough from the Fermi surface. In that 
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